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LONG-TERM  GOALS 

Our  long-tenn  goals  are  to  understand  the  dynamics  of  atmospheric  motions  on  scales  of  order  10  m- 
10  km  in  sufficient  detail  to  be  able  to  provide  a  consistent  subgrid  scale  model  that  will  represent  the 
influence  of  turbulent  transport  across  a  range  of  scales. 

OBJECTIVES 

The  chief  objective  of  the  present  grant  is  to  better  understand  the  physical  processes,  which  control 
cloud  and  circulation  structure  in  the  atmospheric  boundary  layer  and  rate  of  entrainment  of  heat  and 
moisture  across  the  capping  inversion.  This  understanding  will  be  used  to  formulate  a  consistent, 
robust  single-column  model  to  represent  these  processes  in  Mesoscale  models  for  a  broad  range  of 
boundary  layer  conditions  and  forcings. 

APPROACH 

This  research  utilizes  high-resolution  turbulent  transport  codes  and  theoretical  understanding 
developed  under  previous  ONR  support.  Our  principal  approach  is  to  employ  large  eddy  simulations 
(LES)  to  conduct  controlled  numerical  studies  of  the  effects  of  different  boundary  layer  forcings  and 
conditions  (initial  temperature  and  moisture  profiles,  surface  heat  and  moisture  fluxes,  cloud-top 
radiation,  wind  shear,  etc.)  on  the  boundary  layer  dynamics,  cloud  structures  and  entrainment  rates  that 
result.  The  simulations  are  motivated  by  and  compared  with  field  observations  when  available.  The 
simulation  results  are  used  to  develop  and  test  theoretical  models  for  the  basic  physical  processes  at 
work,  in  order  that  these  effects  can  eventually  be  consistently  incorporated  into  lower  resolution 
models. 

WORK  COMPLETED 

During  the  past  year  we  have  concentrated  on  developing  and  testing  buoyancy  flux  models  for  partly 
cloudy  layers  (Lewellen  and  Lewellen,  2003,  LL03  hereafter).  We  developed  an  approximate 
expression  for  buoyancy  flux  formulated  solely  in  terms  of  the  liquid  potential  temperature  and  total 
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water  profiles  and  their  respective  flux  profiles.  The  model  predictions  compare  favorably  with  the 
results  of  an  extensive  set  of  large  eddy  simulations,  including  simulations  of  stratocumulus,  shallow 
cumulus,  and  transitional  behavior  in  between.  This  formulation  was  combined  with  previous  results 
(Lewellen  and  Lewellen,  2002)  on  the  relation  between  cloud-top  entrainment  rate  and  circulation 
structure  to  predict  the  behavior  of  quasi-steady  cumulus-coupled  boundary  layers  as  a  function  of  a 
basic  set  of  physical  input  parameters.  These  predictions  also  compare  quite  favorably  with  LES 
results. 

We  have  also  conducted  a  series  of  simulations  aimed  at  understanding  the  observations  of  nocturnal 
marine  stratocumulus  taken  in  July  2001  during  DYCOMS-II  (the  second  study  of  the  Dynamics  and 
Chemistry  of  Marine  Stratocumulus)  (Stevens  et.  ah,  2003).  Some  of  these  simulations  are  for  a  model 
intercomparison  coordinated  by  the  GCSS  (GEWEX  Cloud  Systems  Studies)  boundary  layer  cloud 
modeling  working  group.  Results  from  a  previous  GCSS  intercomparison,  simulations  of  the  diurnal 
cycle  of  marine  stratocumulus,  were  written  up  and  submitted  for  publication  (Duynkerke  et  al.  2003). 

We  have  participated  in  two  additional  intercomparison  collaborations:  submitting  results  (LES  and 
ID)  of  an  artic  stable  boundary  layer  case  for  the  first  GABLS  (GEWEX  Atmospheric  Boundary  Layer 
Study)  (http://www.gewex.org/gabls.htm)  intercomparison;  and  simulations  of  convective  boundary 
layers  in  a  study  of  the  effects  of  shear  on  entrainment  organized  by  Dr.  E.  Fedorovich. 

RESULTS 

In  a  partly  cloudy  atmospheric  boundary  layer  we  can  formally  write  the  horizontal  mean  buoyancy 
flux  as, 


<w'B’>  =  (1-R)D  +  RW  (1) 

where  D(z)  and  W(z)  represent  the  "dry"  and  "wet"  buoyancy  fluxes  that  would  result  if  the  layer  were 
all  clear  or  all  cloudy,  respectively,  at  height  z,  and  R(z)  represents  the  ratio  of  the  actual  liquid  water 
flux  to  the  liquid  water  flux  that  would  arise  for  completely  cloudy  conditions.  To  good  approximation 
D  and  W  can  be  expressed  as  linear  combinations  of  the  fluxes  of  two  conserved  quantities,  liquid 
potential  temperature  (0i)  and  total  water  (qt).  If  the  vertical  velocity  and  cloud  water  were  completely 
uncorrelated  then  (1)  would  hold  with  R(z)  equal  to  the  mean  partial  cloud  fraction.  This  is  generally 
not  the  case,  however,  particularly  for  shallow  cumulus,  where  the  actual  R  can  exceed  the  cloud 
fraction  by  a  large  factor.  For  this  reason  buoyancy  flux  modeling  in  partly  cloudy  boundary  layers 
has  been  a  longstanding  challenge. 

Motivated  by  LES  results  and  mass  flux  models  of  cumulus  plumes,  we  have  developed  an 
approximate  formulation  of  R  that  can  be  written  in  the  form  (LL03), 

R  ~  (  qisc  -  qib  )/(ssc  -  sb)  (2) 

where  qisc(z)  is  the  liquid  water  of  a  parcel  at  height  z  that  has  near-surface  values  of  9i  and  qt  (with  the 
fluctuations  in  those  values  approximately  accounted  for),  sSL  is  the  difference  between  the  total  water 
and  saturation  mixing  ratios  of  such  parcels,  and  qih  and  sb  are  the  analogous  quantities  for  parcels  with 
the  mean  9i  and  qt  at  the  height  z. 
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The  only  information  we  use  to  evaluate  (2)  are  the  mean  9i  and  qt  profiles.  This  formulation  has  the 
advantage  of  not  requiring  knowledge  of  less  accessible  quantities  such  as  the  mean  liquid  water,  cloud 
fraction,  saturation  variance,  or  vertical  velocity  skewness,  as  in  related  formulations.  Comparison 
with  a  broad  range  of  LES  results  has  shown  that  the  formulation  does  a  good  job  representing  the 
buoyancy  flux,  bridging  between  the  limits  of  layers  that  are  well  characterized  by  joint-Gaussian  pdfs 
(such  as  stratocumulus)  and  cumulus-coupled  layers  with  strongly  correlated  clouds  and  updrafts  and 
very  low  cloud  fractions.  The  formulation  perfonns  particularly  well  for  transitional  layers  in  between 
uniform  stratocumulus  and  well  developed  shallow  cumulus.  The  behavior  in  such  cases  (e.g., 
cumulus  under  stratocumulus)  differs  from  that  predicted  by  the  formulation  of  Bechtold  and  Siebesma 
1998,  which  was  empirically  fit  to  the  limits  of  very  high  or  very  low  cloud  fraction.  In  the  present 
formulation  there  is  not  a  one-to-one  correspondence  between  the  effective  cloud  fraction  R  and  the 
mean  cloud  fraction.  The  underlying  dynamical  regime  is  important  and  is  reflected  in  the  structure  of 
the  9i  and  qt  profiles.  This  conclusion  is  supported  by  the  LES  results. 

This  formulation  for  the  buoyancy  flux  allowed  us  to  extend  our  previous  work  on  entrainment  and 
quasi-steady  structure  in  cloudy  boundary  layers:  first  using  it  to  predict  entrainment  rates  for  partly 
cloudy  cumulus  layers  that  are  in  good  agreement  with  LES  results;  and  second,  using  this  result  to 
predict  the  equilibrium  structure  in  those  layers.  In  Lewellen  and  Lewellen  2002  it  was  argued  that  the 
feedback  between  large-scale  circulation  structure  and  entrainment  into  the  boundary  layer  from  above 
is  a  critical  ingredient  in  detennining  the  quasi-steady  properties  a  layer  will  equilibrate  to,  including 
the  0i  and  qt  profile  shapes.  By  postulating  a  one-parameter  family  of  profiles  for  these  conserved 
variables,  it  is  possible  to  use  the  entrainment  constraints  to  choose  from  among  this  family.  The 
predictions  require  only  a  minimal  set  of  input  parameters  to  be  given  (the  surface  fluxes  of  9i  and  qt, 
their  mean  subcloud  values,  their  respective  jumps  across  the  top  of  the  layer  at  z;,  and  specified 
vertical  profiles  of  any  sources  they  may  have  within  the  layer)  and  the  results  are  in  good  agreement 
with  those  from  LES  for  an  extensive  set  of  cumulus-coupled  cases.  Figure  1  illustrates  examples  of 
the  comparison.  Given  the  limitations  of  the  one  parameter  family  of  9i  profiles  chosen  (which  were 
simply  assumed  linear  within  the  cloud  layer),  the  predictions  for  all  the  profiles  capture  the  basic 
structure  of  the  quasi-steady  LES  results  quite  well.  Much  is  predicted  from  a  minimal  set  of  physical 
inputs,  without  invoking  a  posteriori  features  of  the  LES  results  such  as  the  buoyancy  flux,  cloud 
fraction,  entrainment  rate,  or  temperature  and  moisture  differences  between  the  cloud  and  subcloud 
layers. 

The  success  of  these  predictions  for  many  different  cases  supports  the  validity  of  the  underlying 
physical  balances  assumed  to  produce  the  quasi-steady  equilibrium  for  shallow  cumulus  layers.  This 
evidence  for  a  large-scale  balance  constraining  the  buoyancy  flux  may  in  part  explain  the  results  from 
LES  studies  (e.g.,  Brown  1999)  showing  insensitivity  to  resolution  and  subgrid  model  of  the  mean  flux 
profiles  for  quasi-steady  cumulus  simulations  even  though  the  structure  of  the  individual  cumulus 
plumes  shows  a  clear  sensitivity. 

Our  LES  studies  of  the  first  research  flight  of  DYCOMS-II  seem  to  support  one  of  the  most  interesting 
findings  of  the  observations:  the  maintenance  of  a  solid  stratocumulus  cloud  deck  and  a  coupled  well- 
mixed  boundary  layer  despite  cloud-top  conditions  that  would  seem  to  foster  buoyancy  reversal,  cloud 
breakup  and  decoupling  of  the  layer.  This  result  does  show  some  sensitivity  to  subgrid  modeling  of 
buoyancy  reversal  at  cloud  top.  An  upcoming  GCSS  model  intercomparison  of  simulations  based  on 
this  case  should  address  how  robust  these  results  are. 
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Figure  1.  Mean  vertical  profiles  from  three  sample  large  eddy  simulations  of  quasi-steady  cumulus- 
coupled  boundary  layers  (red,  black,  and  green  lines)  showing:  liquid  water  potential  temperature 
flux,  total  water  flux,  liquid  water  potential  temperature,  total  water,  R  described  in  the  text,  and 
buoyancy  flux.  The  solid  lines  represent  a  horizontal  time  average  of  the  LES  data  and  the  dashed 
lines  the  predictions  discussed  in  the  text  following  from  the  entrainment  constraints  and  given 
values  for  a  limited  set  of physical  inputs.  The  predictions  provide  a  good  representation  of  the  LES 

data  for  all  profiles. 

Results  from  a  previous  GCSS  intercomparison  (Duynkerke  et  al.  2003)  demonstrated  that  LES  does  a 
good  job  reproducing  the  basic  features  of  the  observed  diurnal  cycle  of  marine  stratocumulus, 
including  the  variations  in  liquid  water  path  and  solar-induced  decoupling  of  the  layer  in  the  afternoon. 
Single  column  models  simulating  the  same  case  showed  significant  variations  and  compared  less 
favorably  with  the  observations,  in  large  part  because  of  problems  correctly  predicting  the  cloud-top 
entrainment  rate. 

IMPACT/APPLICATION 

A  valid  model  of  the  marine  atmospheric  boundary  layer  is  needed  not  only  because  its  dynamics  are 
critical  in  transporting  the  surface  fluxes  of  heat  and  moisture  responsible  for  driving  most  atmospheric 
dynamics,  but  also  because  it  determines  the  immediate  environment  within  which  many  Navy 
operations  are  perfonned.  Boundary-layer  clouds  directly  influence  both  the  local  environment  (e.g., 
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for  visual,  radar,  and  laser  communication)  and  the  global  environment  (through  their  impact  on  the 
radiation  balance).  Accurate  representation  of  the  buoyancy  flux,  entrainment  rate,  and 
cloud/circulation  structure  are  critical  ingredients  in  correctly  modeling  the  behavior  of  cloudy 
boundary  layers.  We  expect  this  effort  to  lead  to  improved  subgrid  parameterization  of  these  processes 
in  models  such  as  the  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  System  (COAMPS)  regional 
model  developed  at  the  Naval  Research  Laboratory  (Hodur,  1997). 

TRANSITIONS 

Stevens  et  al.  2003  found  the  entrainment  parameterization  based  on  Lewellen  and  Lewellen  1998  to 
be  one  of  the  few  parameterizations  they  tested  that  was  consistent  with  their  entraimnent 
measurements  taken  during  DYCOMS-II. 

Intercomparison  results  from  the  GCSS  boundary  layer  cloud  modeling  working  group 
(http://www.atmos.washington.edu/~breth/GCSS/GCSS.html)  have  become  a  standard  test  set  for 
development  and  evaluation  of  parameterizations  used  in  numerical  weather  prediction  and  climate 
models. 

RELATED  PROJECTS 

The  LES  code  developed  under  ONR  support  has  been  modified  and  used  to  model  aircraft 
wakes/contrails  for  NASA  (e.g.,  Lewellen  and  Lewellen,  2001),  and  to  model  the  turbulent  interaction 
of  a  tornado  with  the  surface  for  NSF  (e.g.,  Xia  et  al.  2003).  The  use  of  essentially  the  same  LES  code 
on  these  separately  supported  efforts  works  to  the  advantage  of  all  three  projects,  particularly  in 
fostering  numerical  improvements  in  the  efficiency  and  accuracy  of  the  code. 
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